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ABSTRACT 

PSR B0355+54 is one of the handfnl of pnlsars that has been observed with 
both Chandra and XMM-Newton. The analysis of the archival data has revealed 
the pnlsar and a ~ 30" compact nebnla snrronnding it. XMM-iVewton also has 
detected a trail which extends ~6' and similar to the compact nebnla is also 
connter-aligned with the proper motion of the pnlsar. The spectrnm of both the 
pulsar and the extended emission are well described by an absorbed power-law 
model. The measured flux corresponds to an efficiency of converting the spin- 
down luminosity into X-rays in the 2-10 keV band of ~0.01% and ~1% for the 
pulsar and the extended emission, respectively. From the XMM-Newton data we 
have detected pulsations at the expected radio period. The energetic and the 
extent of the extended emission can be explained by a Bow-shock formed by the 
motion of the pulsar through the interstellar medium. 

Subject headings: pulsars:individual (PSR B0355-I-54) — stars:neutron — X- 
rays:stars 


1. Introduction 

PSR B0355-I-54 is a middle age (r ~0.56 Myr) pulsar with a spin-down energy loss rate 
oi E = 4.5 X 10^^ erg s“^, and a relatively close distance of d = l.Od^Q^^g kpc (Chatterjee 
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et al. 2004). It was observed with ROSAT for a total of 19 ks and was detected with a 
non-thermal luminosity (2-10 keV) of 1.2 x 10^^ erg s“^ (Slane 1994). Taking into account 
the extended emission discussed in this paper, one should consider this as only an upper- 
limit to the pulsar’s X-ray emission, considering ROSAT does not have a sufficient angular 
resolution to resolve these features. 

Pulsar wind nebulae (PWNe) are useful probes of pulsars and their surroundings. In 
the case of older (r > 10^ years) isolated neutron stars such as Geminga the PWN is thought 
to be the Bow-shock created by the pulsar proper motion through the interstellar medium 
(ISM) (De Luca et al. 2005; Caraveo et al. 2003). The distance, the energy loss rate, 
and the proper motion of some of the pulsars is well known through radio observations. 
The knowledge of such parameters make pulsar Bow-shocks particularly appealing to study. 
In this letter we report the discovery of an extended emission associated with the pulsar 
B0355-I-54, and investigate the possible implications. 


2. Observations 

2 . 1 . Chandra 

PSR B0355-I-54 was observed using the Advanced CCD Imaging Spectrometer (ACIS) 
instrument aboard the Chandra X-ray Observatory on 2004 July 16. Data were collected 
in the nominal timing mode, with 3.241 s exposures between CCD readouts, and ’’FAINT” 
spectral mode. The standard Chandra screening criteria produced a total usable exposure 
time of 67.13 ks. The ACIS image reveals a point source near the pulsar radio position and 
a relatively faint diffuse emission that extends ~ 30" towards the south-east side, which we 
label as the compact nebula (CN). The detected point source is located at q;=03^58™53!72, 
S= -|-54°13'1379 (J2000.0) which is only 0'.'2 away from the proper motion corrected pulsar 
radio position. This source is detected with an ACIS-S countrate of 0.0031±0.0002 cps in 
the 0.2-10 keV band. The errors quoted in this work are in the la conhdence range, unless 
otherwise noted. 

After removing the contributions from the point sources, we smoothed the image to 
examine the extended emission. The resulting image is given in Figure 1. As can be seen, 
there is excess emission towards south-east of the pulsar. We compared the surface brightness 
of this extended emission to that of the background’s and calculated the countrate of the 
extended emission from an annulus situated at the pulsar position that extends from 3"to 
35". The background countrate was measured from an annulus that extends from 37" to 
70". This annulus was chosen to be co-centric with the extraction region used for the 
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extended emission. The observed extended emission gives an excess of 714±69 counts over 
the background. This corresponds to a net surface brightness of 670±61 counts arcmin”^. 
We also compared the background used for this analysis to a region without any sources 
taken away from the pulsar, towards the chip boundary. This resulted in no discrepancies 
between the surface brightness for the two regions. 


2.2. XMM-ATewton 

PSR B0355+54 was observed with XMM-Newton on 2002 February 10 for 29 ks and 
28 ks with EPIC-MOS and EPIC-pn instruments, respectively. The EPIC-pn EPIC-MOSl, 
and EPIC-MOS2 instruments were operated in the small-window, full-window (imaging), and 
timing modes which provide a temporal resolution of 6 ms, 2.6 s, and 1.5 ms, respectively. 
The medium hlters were used for MOS, the thin hlter was used on pn. The EPIC-pn and 
-MOSl cameras detect a point a source consistent with the Chandra and the radio position. 

The CN detected with Chandra is also visible in the XMM-ATewton data (see Figure 1). 
Although XMM-Xewton data has resolved the extended emission it is difficult to perform 
any spectral analysis on the CN due to relatively low spatial resolution. Consequently, this 
region in EPIC-MOS data (similar for EPIC-pn) contains signihcant amount of counts from 
the pulsar. For example, in an annular region that is centered on the source of interest and 
extends from 10"-35" contains ~35% of the pulsar counts. 

We should also note that in the XMM-Newton data only there is extended emission 
that extends up to ~ 7' that is in the opposite direction of the proper motion, we label this 
feature as the trail (see Figure 1). This structure is aligned with the CN and at an initial 
glance seems fainter. To quantify this, we extracted counts from an elliptical region that 
is aligned with the pulsar’s proper motion and is 3T wide and 5^8 long which was located 
P5 south-east of the pulsar. This region contained 316±48 background-subtracted source 
counts corresponding to a surface brightness of 22.7 counts arcmin”^. When we analyzed the 
Chandra data no such emission was detected. This is not surprising since Chandra is less 
sensitive with respect to XMM-Newton. For example, the Chandra observation of Geminga 
pulsar was unable to conhrm the two tails detected with the XMM-Newton observatory (De 
Luca et al. 2005; Pavlov et al. 2005). 
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3. Timing 

PSR B0355+54 has a period of 0.156 s. When Chandra ACIS is operated in imaging 
mode its temporal resolution is 3.2 seconds which is not sufficient to resolve the periodicity of 
the pulsar. When operated in the small window mode XMM-Newton EPIC-pn camera has a 
timing resolution of 6 ms. This becomes 1.5 ms for the EPIC-MOS2 which was operated in 
the timing mode, where imaging is made only in one dimension to obtain a higher temporal 
resolution. 

In order to check for any pulsations we first extracted source counts from a circle centered 
on the pulsar with a radius 20". Such a small radius was used to reduce the background 
due to the unpulsed extended emission surrounding the pulsar. The extracted event times 
were barycenter corrected. We then calculated the value of the Zf (Buccheri et al. 1983) 
at the predicted pulsar frequency (see Table 1). The obtained value of 7.9 corresponds to a 
probability of chance occurrence of 0.019. We also applied the same criteria to the EPIC- 
MOS2 data. Since the majority of the data in the vicinity of the pulsar is below 0.2 keV we 
took all photons with energies below 2 keV. The Zl for this set at the expected frequency 
was 17.5, which has a probability of 0.00016 of occurring by chance. 

At the radio period we constructed the pulse profiles for both MOS2 and pn event 
times (Figure 2). Both prohles show a double peak nature, however the smaller peak is not 
signihcant, specihcally in the MOS2 data. As was discussed by Becker et al. (2005a) this 
kind of a pulse prohle would not be observed from a neutron star emitting from its hot polar 
caps. Together with spectral data this is suggestive of magnetosphere being the main source 
of X-ray emission for PSR B0355-I-54. 


4. Spectral analysis 
4.1. Pulsar 

For the Chandra observation the pulsar’s energy spectrum was extracted from a 2" 
circular region centered on the pulsar position. Whereas the background was extracted from 
an annular region extending from 37" to 70". The reason for selecting the background away 
from the pulsar is due to the extended emission surrounding the immediate vicinity of the 
pulsar. The source region contains >95% of the source counts. 

The spectrum extracted from the MOSl data was done by selecting all events detected 
in a circle of radius 20" centered on the pulsar position. Using the XMM-Aewton EPIC-MOS 
model point spread function (PSF), 75% of all source counts are within this region. The 
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reason for selecting such a small source region, and in turn fractional encircled energy, was 
to minimize the contamination due to the events originating from the extended emission. 
The background spectrum for the MOSl data was extracted from an annular region centered 
on the pulsar position which extends from 55"-200". For the EPIC-pn data we also used an 
extraction radius of 20" centered on the pulsar. This selection radius includes 75% of the 
point-source flux. Since the point source was close to the edge of the CCD, we extracted 
the background spectrum from a source free region around the pulsar from a circular region 
with a radius of 45". The extracted spectra for PSR B0355-I-54 was binned so as to have 15, 
25 and 20 counts per bin, for ACTS, EPIC-pn and EPIC-MOS, respectively. 

Among the single-component spectral models htted to Chandra spectra a power-law 
model was found to give the best (x^=7.6 for 9 degrees of freedom [dof]) representation of 
the observed energy spectrum of the pulsar. A single blackbody model for the pulsar gave 
neither a statistically nor a physically acceptable £t (x^=21.9 for 9 dof). The parameters 
for the best fit power-law model are given in Table 2. The value obtained for the absorbing 
column is consistent with the value derived through the basic assumption that there is 10 
neutral hydrogen atoms for each free electron, giving a value of Atj^ ~ 2 X 10^^ cm 
The measured X-ray luminosity corresponds to an efficiency of converting the spin-down 
luminosity into X-rays in the 2-10 keV band of 0.009%. 

The EPIC-MOS and -pn spectra were fit simultaneously to get the best fit parameters. 
As was done for Chandra spectrum we tried the single-component models, power-law and 
blackbody. It was found that absorbed power-law model gave the best representation of the 
data (see Table 2). Single blackbody model did not give a statistically acceptable (x^=51.6 
for 33 dof) spectral fit. We also fitted an absorbed power-law model all pulsar spectra (ACIS, 
EPIC-MOS and -pn) simultaneously. This fit yielded a higher photon index but the best 
values for the parameters are consistent with single Chandra and simultaneous XMM-Xewton 
fits (see Table 2). 


4.2. Extended emission 

The spectrum for the CN was selected from an annular region co-centric with the source 
region that has inner and outer radii 3" and 35", respectively. The background for this 
spectrum was chosen to be the same as the pulsar’s spectrum. For the analysis of the trail 
we extracted the spectrum from the region described in § 2.2. The background for this 
spectrum was extracted from an off-source region south-west of the trail. 

The extracted spectra for the CN and the trail were binned so as to have 25 and 50 
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counts per bin, respectively. We limited the generated spectra to the energy interval 0.2-10 
keV. The absorbed power-law model gave a statistically acceptable £t with x^=24.7 for 42 
dof for the CN and x^=15.9 for 22 dof for the trail. Table 2 lists the best fit parameter 
values for the extended emission. The measured values given in Table 2 for the luminosities 
of both the CN and the tail account for only 0.74% of the pulsar spin-down. We also tried a 
thermal bremsstrahlung model to fit the trail. This model also gave a statistically acceptable 
fit with a x^=16.1 for 25 dof. The characteristic temperature of the region is kT = 6 . 5 ^ 3 * 5 ^ 
keV. The absorbing column is consistent with that of the power-law fit, iVpj = 4. 4^34 x 10^^ 
cm“^. The resulting emission measure (f UeUndV) is ~ 1.9 x 10®"^ cm“^, at a distance of 
1.04 kpc. 


5. Discussion 

The Chandra and XMM-iVewton observations of PSR B0355-I-54 have proven the pres¬ 
ence of extended emission associated with the pulsar. This emission has two components 
which are ~30" and ~ 6 ' long (see § 2). The recently detected proper motion direction (see 
Table 1) being counter aligned with the trail together with the observed spectrum being 
well fitted with a power-law model suggests that the observed extended emission may be 
the result of synchrotron radiation from a wind-driven nebula supported by the ram pres¬ 
sure created by the proper motion through the ISM. This hypothesis have been found to be 
consistent with the size and energetics of the X-ray trails of Geminga (Caraveo et al. 2003) 
and PSR B1929-I-10 (Becker et al. 2005a). 

There are two relevant time scales that need to be considered. The first is the Tfiow, 
the timescale for the passage of the pulsar over the length of its X-ray trail. The other one 
is Tsyn, the synchrotron lifetime of the radiation. This is given by the equation (see e.g. 
Gaensler et al. 2002) 

Tsyn = kyr, 

where B^g is the magnetic held in the emission region in units of frG. The first one of these 
can be estimated through the relation Tfiow = "^t/vp, where rt is the extent of the trail which 
is measured to be ~ 6 '. This at a distance of 1.04 kpc corresponds to a linear distance of 
1.8 pc. The proper motion, Vp, is measured by Chatterjee et al. (2004) and is 6llg^ km 
s“^. This, together with the linear extent of the trail, gives the time it took PSR B0355-I-54 
to travel across the trail as ~ 34 kyr. Assuming these two time scales are comparable the 
estimate for the magnetic held around the emitting region is B^g for an X-ray photon of 
energy hu = 1 keV. We should note that this calculation ignores the possible inclination of 
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the pulsar’s proper motion with respect to the line of sight. 

Assuming that the ISM in the vicinity of the pulsar carries a magnetic held that does 
not differ signihcantly from that of the average of our galaxy which is ~ 4 /iG (Beck et al. 
2003). This value should be the characteristic magnetic held for the trail where the X-ray 
emitting particles are no longer in the shocked region. For the CN, where the electrons are in 
the shocked region, we need to consider the change in density due to the shock and hence the 
change in the magnetic held. We can use standard hydrodynamical arguments to show that 
the density in the shock region is given by the relation Pgjiock ~ 4PlSM> which would imply 
a magnetic held of i^ghock ~ -^shock ~ given the average ISM magnetic 

held. This value, given the large uncertainties, is in agreement with the value found here. 

The termination shock radius is given by the balance of the ram pressure between the 
wind particles and the ISM, i.e. 

E 1 2 

where p is the density of the ISM and Vp is the proper motion velocity. From here one can 
estimate the radius of the termination shock as ~ 3 x cm, where n 

and are the number density in the ISM in units of cm“^ and the pulsar space velocity 
in units of 100 km s“^, respectively. Adopting n = 1 cm“^ the shock termination radius 
becomes 1.1x10^^ cm. This corresponds to an angular size of ~7" at a distance of 1.04 kpc. 
This value is consistent with the data given the observational constraints. 

In order to check the energy dependence of the extent of the CN we constructed two 
background-subtracted linear prohles along the proper motion direction, which contain pho¬ 
tons from the intervals 0.2-2 keV and 2-10 keV (Figure 3). For comparison also the point 
source contribution is plotted. Although the structure of the CN is not identical in both 
energy ranges the overall extent within the statistical uncertainties are equal. This is consis¬ 
tent with the identihcation of this region as the termination shock. This region is expected 
to show a constant X-ray spectrum across its extent. Due to the pressure difference between 
the regions ahead and behind the pulsar’s motion the termination shock is not of uniform 
radius around the pulsar and is extended along the proper motion direction (Bucciantini 
et al. 2005). The larger structure seen by XMM-Newton (trail) can be thought of as the 
emission originating from the wind particles that flow around the edge of the shocked region. 
The magnetic held in this region should be an order of magnitude less than the nebular mag¬ 
netic held in the shocked region, if equipartition holds. The two components of the extended 
emission could be explained by the Bow-shock theory. 

PSR B0355-I-54 is yet another pulsar with known proper motion and a collimated outhow 
in the opposite direction. Such how, also observed in Geminga, could be interpreted as a 



pulsar jet (Pavlov et al. 2005). This is reminiscent of the argument proposed for the Vela 
pulsar (Markwardt & Ogelman 1995) where the outflow in opposite direction is proposed as 
the mechanism to boost pulsar velocities in contrast to birth-kicks. 
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Table 1. Properties of PSR B0355+54 derived from radio observations 


Parameter Value 


Frequency (Hz). 6.39458383 

Frequency derivative (Hz s“^). -1.797906x10“^^ 

Epoch (MJD). 52315.6 

Spin-down age (10® yr). 0.56 

Spin-down energy (10®"^ ergs s“^). 4.5 

Inferred Magnetic Field (10^^ G). 8.4 

Dispersion Measure (pc cm“^). 57.2 

Distance (kpc). 1.04 

Position (J2000) 

a. 03^58'^53!72 

S . +54°13'13"9 

Proper motion (mas yr“^) 

/i„ . 9.20±0.18 

fxs . 8.17±0.39 


Note. — Distance, position, and proper motion are all taken 
from Chatterjee et al. (2004). 
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Fig. 1.— Gray-scale plots of the 0.2-10 keV image centered on the pulsar PSR B0355-I-54 
(left) by XMM-Newton EPIC-MOS and (right) Chandra ACIS. The X-ray point source 
contribution is removed, the images are smoothed with a Gaussian of width 2". Also shown 
are the isophoto contours. On the left panel the box shows the zoomed in region shown on 
the right panel (Chandra image). The direction of the pulsar proper motion is shown on the 
right panel with a black arrow and the pulsar position is indicated by a white cross (x). The 
length of the proper motion vector is the distance the pulsar would travel in 1000 years. 










Counts/bin Counts/bin 
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Phase 

Fig. 2.— X-ray pulse profiles of PSR B0355-I-54, constructed from EPIC-MOS2 top and 
EPIC-pn bottom. The pulse fractions are given in the hgure. Notice the double-peak nature 
of the pulse prohles. The peaks are ~0.3 of a period apart. Two cycles are shown or clarity. 



Table 2. 


Spectral fit results for PSR B0355+54. 


Comment 

xV(dof) 

P 

% 

1 

O 

7-2-10 

-^x 




lO^^cm ^ 

10 ’^"’erg cm ^ s ^ 

lOSierg s“’ 

Chandra 

Pulsar 

7.6/9 

1 9A+^-^^ 
-‘-•^^-0.24 

o.7ol^:?« 

q 1 1 +0.40 
'J.-‘--‘--0.43 

0.40 

CN 

24.7/42 

1 44+0-26 
-‘-•++-0.34 

6 .1011:81 

12.311:3'’ 

1.61 

XMM-JVewtou’’ 

pn-|-MOS 

37.3/32 


4 00+2-62 

11-3111J 

1.47 

pn+MOS+ACIS 

50.7/43 

i.4i«:J° 

q qo+l.OO 

O.OZr_^ 

3.i4lS:g 

0.41 

Trail 

15.9/22 

l-844;5? 


13.51^:^ 

1.72 


Note. — The quoted errors are la 

“^The value quoted is the integrated flux in the annular region where the spectra for 
the extended emission was extracted. 

’’The hrst and second row tabulate the best £t parameters for the spectral fit performed 
simultaneously on data extracted from EPIC-pn and -MOS and EPIC-pn, -MOS, and 
ACTS, respectively. The third row shows the values for the absorbed power-law best fit 
parameters. 
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Fig. 3.— The background-subtracted linear profile of the CN, along the proper motion 
direction as seen by Chandra in the energy range 0.2-2 keV {filled circles) and 2-10 keV 
{diamonds). The solid line is the contribution from the pulsar B0355-I-54 in the 0.2-2 keV 
range. 



















